function []=chapter3_1run()
    t = [0.2,0.4,0.6,0.8,1];
    h = [0.1,0.01,0.001];
    for i = 1:length(h)
        for j = 1:length(t)
            %LaxWendrof(t(j),h(i));
            
        end
    end
    for i = 1:length(h)
        for j = 1:length(t)
            %conserve(t(j),h(i));
            t(j),h(i)
        end
    end
    for i = 1:length(h)
        for j = 1:length(t)
            nonconserve(t(j),h(i));
            t(j),h(i)
        end
    end
    
            
end